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Abstract 

The influence of the boundary conditions used in the Wigner-Seitz 
approximation applied to the neutron star inner crust is examined. 
The generalized energy functional method which includes neutron and 
proton pairing correlations is used. Predictions of two versions of the 
boundary conditions are compared with each other. The uncertainties 
in the equilibrium configuration (Z,Rc) of the crust, where Z is the 
proton charge and Rc the radius of the Wigner-Seitz cell, correspond 
to variation of Z by 2 - 6 units and of i?c) by 1 - 2 fm. The effect 
of the boundary conditions is enhanced at increasing density These 
uncertainties are smaller than the variation of Z and Rc coming from 
the inclusion of pairing. The value of the pairing gap itself, especially 
at high density can depend on the boundary condition used. 

PACS : 26.60.+c,97.60.Jd,21.65.+f,21.60.-n,21.30.Fe 

In the last two decades the interest on the structure of the neutron star 
inner crust has been stimulated by the increasing number of observational 
data on the pulsar glitches. The latter are commonly explained in terms of 
the dynamics of superfluid vortices within the inner crust of neutron stars 
(see [Tj and Refs. therein). By "inner crust" one usually indicates the part 
of the shell of a neutron star with sub-nuclear densities O.OOlpo < P ^ 0.5po, 
where po is the normal nuclear density. According to present-day ideas, the 
bulk of the inner crust consists mainly of spherically symmetrical nuclear- 
like clusters which form a crystal matrix immersed in a sea of neutrons and 
virtually uniform sea of electrons. Such a picture was first justified micro- 
scopically in the classical paper by Negele and Vautherin (NV) j2] within the 
Wigner-Seitz (WS) approximation. Up to now, the WS method remains to 
be quite popular in this field. Only recently a more consistent band theory 
was developed for the deep (high density) layers [H], where the "lasagna" or 
"spaguetti" structure of the crust matter is supposed to be favored, and for 
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the outer (low density) layers as well As far as the band theory is quite 
complicated, the WS method is usually considered as the most practical one 
for systematic investigation of the inner crust structure in the whole density 
interval. 

NV used for describing the matter of a neutron star crust a version of 
the energy functional method with density dependent effective mass m*(p). 
In fact, it is very close to the Hartree-Fock method with effective Skyrme 
forces. For a fixed average nuclear density p, the nuclear (plus electron) 
energy functional is minimized for the spherical WS cell of the radius Re- 
A cell contains Z protons (and electrons) and N = A — Z neutrons {A = 
(47r/3)-R|?p). In addition, the /3-stability condition. 



Atn - (/ip + fJ-e) = 0, (1) 

has to be fulfilled, where pp and pe are the chemical potentials of neutrons, 
protons and electrons, respectively. The minimization procedure is carried 
out for different values of Z and Re- The equilibrium configuration (Z, Rc) 
at the considered density corresponds to the absolute minimum in energy 
among all these possible configurations. 

Application of the variational principle to the NV energy functional for 
a WS cell results in the set of the Shrodinger-type equations for the single 
particle neutron functions 0a (r) = Rnij{r)^ijm{'n), with the standard nota- 
tion. The radial functions Rnij{r) obey the boundary condition (BC) at the 
point r = R(.. There exist different kinds of the BC. NV used the following 
one: 



RniAr = Rc) = (2) 

for odd /, and 



dR.. 



nlj 



dr 



0, (3) 



for even ones. Let us denote it as BCl. The use of this BC has been partly 
justified by physical considerations in NV, but the dependence of the results 
on the BC has never been discussed in detail. It is the purpose of the paper to 
study this problem at a quantitative level and to establish the corresponding 
uncertainity, which is inherent to the WS method applied to neutron star 
crust. For this aim, we compare results obtained for the BCl with those 
found for an alternative kind of the BC (BC2) when Eq. ^ is valid for even 
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/ whereas Eq. (jH)), for odd ones. In principle, two additional kinds of the BC 
exist when Eq. or Eq. is used for any /. As it was noted by NV, these 
BC have an obvious drawback for the case of the neutron star inner crust as 
far as they lead to an unphysical irregular behavior of the neutron density 
p„(r) in vicinity of the point r = Rc. Indeed, Pni'^) vanishes in this point in 
the first case and has a maximum in the second one. On the contrary, Pn{r) 
is almost constant nearby the point r = in the case of the BCl or BC2. 

It should be noted that the pairing effects were not taken into account 
in [2| since it was supposed that they are not important for the structure of 
the crust. The reason of such an assumption is the rather small contribu- 
tion of the pairing effects to the total binding energy of the system under 
consideration. Recently, we have generalized the NV approach to describe 
the inner crust by explicitly including the neutron and proton pairing corre- 
lations ini 13 El in a self-consistent way. It turned out that in the whole 
interval of p the equilibrium configuration {Z, Rc) changes significantly due 
to pairing. 

We used the generalized energy functional method jS] which incorporates 
the pairing effects into the original Kohn-Sham (KS) ^0] method. In this 
approach, the interaction part of the generalized energy functional (GEF) 
depends, on equal footing, on the normal densities Pn, Pp, and the abnormal 
ones, Un, Up, as well: 



where Sint is the GEF density. It is the sum of two components, the normal 
and the anomalous (superfiuid) ones: 



where t = n,p is the isotopic index. Just as in the KS method, the prescrip- 
tion m* = m holds to be true. 

To describe the central part of a WS cell with the nuclear cluster inside 
we used the phenomenological nuclear GEF E^^ by Fayans et al. which 
describes properties of the terrestrial atomic nuclei with high accuracy. For 
describing neutron matter surrounding the cluster we used a microscopic 
energy functional for neutron matter based on the Argonne NN potential 
vl8 [TT]. The ansatz of |71IH] for the complete energy functional is a smooth 
matching of the phenomenological and the microscopic functionals at the 
cluster surface: 



£ip^{v), u^{v)) = £'\pAr), z..(r))F^(r) +f ™(p.(r), z/.(r))(l - F„(r)), (6) 




(4) 



'int 



■normyPr) £an\pT^ 



(5) 
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where the matching function Fm{r) is a two-parameter Fermi function: 

F^r) = (1 + exp((r - Rj/d^)y\ (7) 

Eq. (jU)) is apphed both to the normal and to the anomalous components 
of the energy functional. After a detailed analysis, the matching parameters 
were chosen as follows. The diffuseness parameter was taken to be equal to 
dm=0.3 fm for any value of the average baryon density of the inner crust and 
for any configuration {Z,Rc). As to the matching radius Rm, it should be 
chosen anew in any new case, in such a way that the equality 

Pp{Rj = O.lppiO) (8) 

holds. In this case, on one hand, for r < R^ neutrons and protons coexist 
inside the nuclear-type cluster, and the use of a realistic phenomenological 
energy functional seems reasonable. On the other hand, at r > Rm one can 
neglect the exponentially decaying proton "tails" and consider the system as 
a pure neutron matter for which an adequate energy functional microscopi- 
cally calculated can be used. The same matching parameters were used for 
normal and anomalous parts of ©. As far as practically all the protons are 
located inside the radius Rm, the matching procedure concerns, in fact, only 
neutrons, protons being described with the pure phenomenological nuclear 
GEF. 

It is worth to mention that for neutron matter region, the ansatz is, 
in fact, the LDA for the microscopic part of the GEF. As it is commonly 
known, the LDA works well only provided the density is smoothly varying, 
whereas it fails in the surface region with a sharp density gradient. The 
above choice of the matching procedure and the values of the parameters 
guarantees that this region of a sharp density variation is mainly governed 
by the phenomenological nuclear part of the GEF which "knows how to deal 
with it" . 

For the microscopic part of the normal component of the total energy 
functional (0) we follow refs. [71 |H] and take the EOS of neutron matter 
calculated in jj^l with the Argonne f 18 potential on the basis of Brueckner 
theory, taking into account a small admixture of 3-body force. Its explicit 
form could be found in the cited articles. The microscopic part of the anoma- 
lous component of the GEF in ^|8j was calculated for the same vl8 potential 
within the BCS approximation. 
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As calculations of El 13 IE] have shown, the pairing correlations influ- 
ence the equilibrium values {Z, Rc) significantly. To explain this effect, it is 
instructive to analyze the /3-stability condition (P). As far as electrons in the 
inner crust of a neutron star are ultra-relativistic, the following relation is 
valid: He ^ (97rZ/4) ^/^/i2c- By substituting it into Eq. (P), one finds 

^-^(/^n-/^p)'^c- (9) 

The influence of pairing on the chemical potentials and Hp is much stronger 
than that on the total binding energy. Their variation may be of the order 
of the gap value A ^1-2 MeV. Such a variation of fin or fip may lead to a 
sizable change of the equilibrium value of Z as far as the difference — fip) 
is raised to the third power in Eq. (jH)). The estimate of the change of Z 
induced by this variation is as follows: 6Z=3Z6{fin — fip)/{fin — f^p)- For 
average values of fcp, the difference /i„ — yUp ~ 50 ^ 70 MeV, hence 5Z could 
reach several units of Z. An additional change of the Z value may appear 
due to a variation of Re- 

Besides, as it can be seen in Fig. 1, the binding energy Eb is rather 
fiat function of Z and different local minima Eb{Z) have often close values 
of Eb- Therefore their relative position may change after switching off the 
pairing since in general the corresponding contribution to is an irregular 
function of Z. Such a situation does often occur within the WS approach, 
especially for high density values, due to the shell-type effect in the single- 
particle neutron spectrum. An example is discussed below. 

In the calculations of [21011113 the NV boundary condition, BCl, was 
used. Here we repeat the analysis for the case of the boundary condition BC2. 
Results for the binding energy are shown in Fig. 1. Just as in [3 E] only 
even values of Z are used. The detailed comparison is made for kp=0.8 fm~^. 
Although the two curves Eb{Z) are quite different, the positions of local 
minima for the BCl and BC2 are close to each other, the distance being equal 
to 2 or 4 units of Z. What is of primary importance, the relative position of 
local minima for BC2 is the same as for BCl. In particular, the positions of 
the absolute minimum almost coincide (Z=52 for BCl and Z=54 for BC2). 
These observations permit us to simplify calculations for other values of kp. 
In the case of the BC2, we limit ourselves mainly with the analysis of a 
vicinity of the absolute minimum for the BCl. The neighborhood of other 
local minima was analyzed only in the case if they have values of Eb{Z) 
close to that corresponding to the absolute minimum. It turned out that 
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Figure 1: Binding energy per nucleon for various in the case of the 
BCl (sohd circles connected with the sohd hnes) and the BC2 (open circles 
connected with the dotted lines). 
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there is no value of kp for which the relative position of a local minimum 
and of the absolute one for the BCl and BC2 is different. In addition to 
systematic calculations for kF=0.6 -r- 1.2 fm~^, we made an extra one for a 
small density, A;f=0.2 fm~^, in vicinity of the neutron drip point. In the last 
case, two curves corresponding to BCl and BC2 practically coincide. For all 
other values of k^ the absolute minima are shifted by 2, 4 or even 6 units of 
Z. 



Table 1: Comparison of properties of equilibrium configurations of the WS 
cell for two different kinds of the boundary condition 



kp 


Z 


/?,. fm 


En. MeV 


//n. MeV 


Ai,-. MeV 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


BCl 


BC2 


0.2 


52 


57.18 


57.10 


-0.9501 


-0.9483 


0.1928 


0.1942 


0.04 


0.05 


0.6 


58 
56 


37.51 
36.97 


37.48 
36.95 


2.1516 
2.1563 


2.1596 
2.1572 


3.2074 
3.2173 


3.2226 
3.2193 


1.92 
1.91 


1.89 
1.89 


0.7 


52 
48 


32.02 
31.16 


32.04 
31.14 


2.7908 
2.7924 


2.7989 
2.7856 


3.9876 
4.0069 


4.0107 
3.9873 


2.30 
2.29 


2.25 
2.32 


0.8 


42 
44 


26.90 
27.29 


26.91 
27.30 


3.4373 
3.4435 


3.4471 
3.4319 


4.8454 
4.8553 


4.8561 
4.8198 


2.56 
2.53 


2.45 
2.56 


0.9 


24 
22 


20.26 
19.87 


20.30 
19.70 


4.1123 
4.1141 


4.1169 
4.1104 


5.7340 
5.7861 


5.7986 
5.7170 


2.64 
2.62 


2.51 
2.54 


1.0 


20 
24 


16.69 
18.29 


16.90 
18.22 


4.8210 
4.8479 


4.8522 
4.8231 


6.8525 
6.8446 


6.7424 
6.8920 


2.02 
2.52 


2.52 
2.29 


1.1 


20 
26 


14.99 
16.75 


15.33 
17.08 


5.5765 
5.6677 


5.6733 
5.6100 


7.4288 
7.9680 


8.0446 
8.5398 


1.32 
2.28 


2.32 
2.02 


1.2 


20 
26 


13.68 
15.21 


13.95 
14.89 


6.4225 
6.6639 


6.6762 
6.4587 


8.5814 
9.0825 


9.1898 
9.3413 


1.21 
1.25 


1.56 
0.86 



Comparison of different properties of the equilibrium configuration of the 
WS cell for various values of kp in the case of the BCl and BC2 is presented 
in Table 1. There are two lines for every value of kp. The first one is given 
for the Z value corresponding to the minimum of Eb in the case of the BCl, 
the second one, for the BC2. The only exception is A;f=0.2 fm~^ when these 
two values of Z coincide. In the last two columns, the average value Ap of 
the diagonal matrix element of the neutron gap at the Fermi surface is given. 
The averaging procedure involves 10 levels above /x„ and 10 levels below. 
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One can see that the influence of the BC is enhanced at increasing values 
of kp. EspeciaUy strong variation of Ap and yU„ takes place in the cases of 
kp—l.l fm~^ and A;f=1.2 fm~^ (Fig- 2). To illustrate the influence of the 




5 10 15 

r, fm 

Figure 2: The neutron gap for A';f=1.1 fm~^, Z=20 and Z—26, in the case 
of the BCl (solid hues) and the BC2 (dashed lines). 

BC to the neutron gap in the first case, the gap function A„(r) is drawn for 
both values of Z and both kinds of the BC. The most strong variation of the 
gap occurs in the case of Z=20. To understand the reason of such strong 
effect, we draw the neutron single particle spectrum ex for this value of Z 
in Fig. 3 for the BCl (the left half of the figure) and the BC2 (the right 
one). The position of the chemical potential /i„ is shown with dots. The two 
spectra are absolutely different. The reason is the shift Aex of each A- level 
going from BCl to BC2. The value of this shift is approximately equal to a 
half of the distance between two neighboring levels with the same (/, j), the 
sign of the shift being opposite for even and odd /. The absolute value of 
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Figure 3: The neutron single-particle spectrum sx for kp—l.! fm ^, Z—20, 
in the case of the BCl (left) and the BC2 (right). 



the shift is proportional to l/R^ and grows at increasing values of kp. The 
corresponding shifts arc shown in Fig. 3 for two states, 2ji3/2 and 1^23/2, 
which are the neighbors of /i„ in the BCl case. On average, the spectrum is 
quite dense, however in both cases there is a shell type structure with rather 
wide intervals between some neighboring levels. If one deals with a big inter- 
level space in vicinity of fin, as in the BCl case in Fig. 3, one usually obtains 
a dense set of levels in this region when going to the opposite kind of the BC. 
In the case of the BC2, big intervals are far from the Fermi surface and do not 
influence signiflcantly the value of the neutron gap. On the contrary, in the 
case of the BCl /x„ is situated just inside such an interval that suppresses the 
gap significantly. In principle, the neutron gap could vanish if the interval 
was wider. 
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As we see, there are internal uncertainties inherent to the WS method 
apphed to the neutron star inner crust which originate from the kind of the 
BC used. Only in the case of very small density nearby the neutron drip point 
predictions of the BCl and BC2 versions are practically identical. If we deal 
with kp > 0.6 fm~^, the uncertainty in the equilibrium value of Z is between 
2 and 6 units, growing with increase of kp. The uncertainty in the value of 
Rr rule, about 1 fm and only for kp=l.l fm ^ it turns out to be about 

2 fm. However, the value of these uncertainties is less than the variation of 
the equilibrium configuration {Z,Rc) connected with the pairing effects [S]. 
In the case of high densities, kp > 1 fm~^, the most important uncertainty 
occurs in the value of the neutron gap A„. It originates from the shell effect 
in the neutron single-particle spectrum which is rather pronounced in the 
case of big kp and, correspondingly, small Rc values. 

The authors thank N.E. Zein for valuable discussions. This research was 
partially supported by the Grant NS-1885.2003.2 of The Russian Ministry 
for Science and Education. 
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